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Abstract The continuous high spatial-resolution Doppler observation of the Sun 
by Solar Dynamics Observatory / Helioseismic and Magnetic Imager allows us to 
compute helioseismic k~uj power-spectrum diagram using only oscillations inside 
a sunspot. Individual modal ridges can be clearly seen with reduced power in the 
k-Lu diagram constructed by use of 40-hour observation of a stable and round 
sunspot. Comparing with the k-uj diagram obtained from a quiet-Sun region, 
inside the sunspot the /-mode ridge gets more power reduction than p-mode 
ridges, especially at high wavenumber. The p-mode ridges all shift toward lower- 
wavenumber (or higher- frequency) areas for a given frequency (or wavenumber) , 
implying an increase of phase velocity beneath the sunspot. This probably results 
from acoustic waves' travel across the inclined magnetic field of the sunspot 
• penumbra. Line-profile asymmetries exhibited in the p-mode ridges are more 

significant in the sunspot than in quiet Sun. Convection inside the sunspot is 
also highly suppressed, and its characteristic spatial scale is substantially larger 
than the typical convection scale of quiet Sun. These observational facts demand 
a better understanding of magnetoconvection and interactions of helioseismic 
waves with magnetic field. 
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1. Introduction 

Analysis of helioseismic power-spectrum diagrams, also known as k-tu diagrams, 
obtained from solar oscillation signals, is a useful way to derive solar atmospheric 
properties, interior structures, and subsurface dynamics (e.g. [Rhodes et aO 119971 
ISchou et ail I1998p . While many global-scale properties of the Sun were inferred 
by analyzing k-ui diagrams obtained from the global-scale observations, some 
local-scale properties were derived by analyzing k-uj diagrams obtained from 
oscillations in local areas, known as ring-diagram analysis ( [Hill, 1998| . Solar 
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areas containing active regions have been studied using ring-diagram techniques 
to infer int erior structures (|Basu, Antia, and Bogart, 2004 Baldner, Bogart, and 



Basu, 2011) and subsurface flow fields ( |Haber et ai, 2004[|Konim et ai, 2005P 
of active regions. However, the k-uj diagrams used in these analyses are largely 
composed of oscillations outside of sunspots rather than those inside sunspots, 
hence represent properties of large areas instead of properties of confined sunspot 
areas. Other local helioseismology techniques, e.g.^ time-distance helioseismology 
dPuvall et al, 1993b|[K"osovichev, Duvah, and Scherrer, 2000l|Gizon et al, 2009[ ), 
also need to apply filters over k-oj diagrams although not to analyze the k-oj 
diagrams directly. 

Despite the scientific interest and importance of constructing k-u diagrams 
using only oscillations inside sunspots, this has not been successfully achieved 
due to c ither poor spatial resolution or limited temporal coverage. Penn and 
LaBonte (|1993p attempted to compute k-uj diagrams of two sunspots using 
observations from Mees Solar Observatory, but were not able to obtain distinct 
modal ridges. Using Hinode Ca H observations, [Nagashima et a/.|(P"007p studied 
acoustic power maps inside a sunspot, but did not analyze its k-uj diagram 
due to the short temporal coverage. Employing Hankel's decomposition, Braun, 
Duvall, and LaBonte ^^^TIT analyzed the k-uj diagrams made from acoustic 
waves traveling into and out from sunspots separately, but these diagrams were 
not computed using oscillations observed directly inside sunspots. SOHO/MDI 
( [Scherrer et al., 1995[ ) full-disk data did not have sufficient spatial resolution 
to construct k-uj diagrams using only signals inside sunspots, while its high- 
resolution data lacked sufficient coverage time. The Helioseismic and Magnetic 
Imager ( [Scherrer et al., 2012[ [Schou et al., 20T2| ) onboard Solar Dynamics Ob- 
servatory (SDO/HMI) provides continuous observations of the Sun with a high 
spatial resolution and steady temporal cadence, and this allows us to perform 
a k-Lu diagram analysis for sunspots. In this article, we report some interesting 
properties found in the k-uj diagram of a sunspot. We introduce our data-analysis 
procedure in Section 2, and present our analysis results in Section 3. We discuss 
these results in Section 4, and give conclusions in Section 5. 



2. Data Analysis 

SDO/HMI has continuous observations of the full solar disk except during some 
short eclipse periods, with a temporal cadence of 45 seconds and a spatial 
resolution of 1.0 arcsec. Dopplergrams observed by HMI, primarily used for 
helioseismology studies as recently demonstrated by IZhao et~ai\ (|2012p , are used 
in this study. The sunspot located inside NOAA AR 11092 is selected for the k- 
io power-spectrum analysis. This sunspot, as displayed in Figure [I] was located 
at the latitude of 12.4°N, and remained stable and round in shape during the 
period of analysis. Magnetic field surrounding this sunspot was relatively simple. 
A period of 40 hours, covering 16:00 UT 2 August through 08:00UT 4 August, 
2010, is selected for analysis, and this is roughly corresponding to the sunspot's 
central meridian crossing, from 8.8° East of the central meridian to 8.8° West of 
it. This 40-hour period is then divided into five eight-hour segments, and each 
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Figure 1. Continuum intensity image (a) and line-of-sight component of the magnetic field 
(b) for AR 11092. The magnetic-field strength is displayed with a range of —2000 to 500 Gauss. 
Oscillations observed inside the black circle in panel (a) are used to calculate the k-ui diagram 
for the sunspot. 



data segment is analyzed separately. The results obtained from each segment are 
then averaged to give the final results. Every eight-hour data segment is tracked 
with the sunspot's rotation speed as given by Howard (1990), and is remapped 
using Postel's projection with the remapping center at the sunspot center and 
with a sampling resolution of 0.03° pixel" ^, slightly ovcrsampling in most areas 
except near the solar-disk center. 

Fourier transform is used to convert data from space - time domain to Fourier 
domain, and k-tu diagram is constructed in the three-dimensional Fourier do- 
main. To avoid possible complications of Fourier analysis due to the irregular 
shape of the sunspot boundary, we choose to use the oscillations inside a circular 
boundary as denoted by the black circle in Figure [TJi. When performing Fourier 
transforms, we keep the entire box size unchanged as shown in Figure [T] and fill 
up the area outside of the circle with zeroes, and this is a normal practice known 
as zero-padding in Fourier analysis. The area near the boundary of the inside 
and outside of the circle is tapered with a cosine bell to avoid sharp transition 
that may give spurious signals. For the k-uj diagram made following the above 
procedure, the power at low-wavenumber areas gets contaminated due to the 
small box size and the circular structure with data surrounded by zeroes. Thus, 
we do not take signals lower than k = 0.3 rad Mm^^ into the analysis that 
follows. 

To better understand the analysis results, it is useful to do similar analyses 
over quiet-Sun regions for comparison. Here, we select a piece of quiet Sun 
that has a same disk location and a same temporal duration as the analyzed 
sunspot covering the period of 16:00 UT 17 August through 08:00UT 19 August, 
2010, fifteen days after the sunspot analysis period. The exactly same analysis 
procedures are applied on these data segments, and results are used as references 
to compare with the results from the sunspot. 
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Figure 2. Panels (a) and (b) display the k-u) diagrams obtained from the sunspot and from 
the quiet-Sun region, respectively. Color scales in these two panels are the same. Panel (c) 
shows the ridge-peak locations, obtained by fitting both k-u) diagrams, to better illustrate the 
relative locations of the modal ridges. 



3. Results 

Figure and b display the k-uj diagrams calculated from the sunspot region 
and the quiet-Sun region, respectively. Individual modal ridges can be clearly 
seen in both diagrams. As expected, the helioseismic power of each ridge is 
much weaker in the diagram obtained from the sunspot than that from the 
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Figure 3. (a) Phase velocities calculated from fitted power ridges pi — P3, with solid lines 
for the quiet-Sun region and dashed lines for the sunspot. (b) Perturbations of phase velocity 
obtained from the sunspot relative to that obtained from the quiet-Sun region. 



quiet-Sun. Compared with the p-mode (acoustic mode) power reductions for 
the sunspot, the power reduction in the /-mode (surface-gravity mode) ridge is 
more substantial, and the length of this ridge seems more shortened due to the 
higher power reduction at higher wavenumber. For both k-to diagrams, we fit the 
Pi —p3 modal ridges by peak- finding to better illustrate how the location of each 
individual ridge obtained from the sunspot differs from that obtained from the 
quiet-Sun region. The fitting results are shown in Figure[5J:. It is difficult to get a 
reliable fitting for the /-mode ridge of the sunspot. The three fitted modal ridges 
obtained from the sunspot all shift to the lower wavenumber (or alternatively 
speaking, higher frequency) side relative to the ridges obtained from the quiet 
Sun for a given frequency (or wavenumber). The shift amount is not uniform 
along each individual ridge, and the largest shift is up to 0.3 mHz. 

Based on the fitted modal ridges shown in Figure HJ;, we calculate phase 
velocity Wph = lo/k along each ridge of both diagrams and display the results in 
Figure [3^. The derived phase velocity in the sunspot region is often larger than 
in the quiet-Sun region by up to 3 km s~^. Figure |3Jd shows the phase- velocity 
perturbation for the sunspot relative to the quiet-Sun region (Jwph/i'ph for each 
modal ridge. The calculated (5?;ph/wph, although large in error bars, mostly fall 
in range of 0.02 to 0.06. 

Comparisons of the ridge shape and location at selected frequencies are dis- 
played in Figure m where the curves are obtained by averaging a 0.2-mHz wide 
band in the /c-w diagrams of Figure [2^ and b. Figure H] again shows that the 
helioseismic power of the /-mode ridge is more suppressed in the sunspot. All 
of the p-mode ridges show clear shifts toward lower wavenumber, and the shifts 
are more significant in higher frequencies. The amount of ridge shift can be 
as large as half of the ridge- width for pi, p2, and pz ridges at the frequency 
of 5.0 and 6.0 mHz, as can be seen in Figure SJ; and d. It can also be found 
in Figure 21 that the modal ridge line-profile asymmetry is more significant in 
the power spectrum from the sunspot than from the quiet-Sun region. It is 
well known that the modal ridges obtained from Doppler observations exhibit 
line-profile asymmetries ( |Duvall et ai, 1993a[ ), and this asymmetry is believed 
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to be caused b y correlated noise (jNigam et al., 1998 Georgobiani, Stein, and 



Nordlund, 2003). It is not immediately clear whether the more prominent line- 
profile asymmetries in sunspots can be explained by more substantial correlated 
noise observed in these areas. On the other hand, recent studies (e.o. lChou et al.\ 
I2009|) found that the power absorption and local suppression of acoustic waves 
inside sunspots increase with wavenumber, and the more prominent line-profile 
asymmetry may possibly be due to this effect. 

Figure [S] shows comparisons of ridge shapes and locations at some selected 
wavenumbers. For each given wavenumber, the pi~pz ridges all display shifts 
toward higher frequency, and the amount of shift can also be up to half of the 
ridge width. The significant line-profile asymmetry is also clear in the p-mode 
ridges obtained from the sunspot region. For the /-mode ridge, it seems that the 
ridge shifts toward lower frequency in panels (c) and (d) of Figure [Sj but that 
may just be due to the enhancement of background noises near these frequencies. 

It is also interesting to compare convection power of the sunspot and the quiet- 
Sun region. The convection powers, primarily below 1.5 mHz, can be seen in the 
lower part of panels (a) and (b) of Figure [H Following procedures employed 
by IChou et al.\ (|1992p and [Georgobiani et al\ (|2007p , we calculate convection 
power as a function of wavenumber by integrating power below the straight line 
separating oscillation power and convection power, i.e. v = l.OSfc, where v is in 
unit of mHz and fc in unit of rad Mm^^. Figure [5] shows the comparison of con- 
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vection powers obtained from both k-uj diagrams. Not surprisingly, convection 
is greatly suppressed inside the sunspot. The location of the convection peak 
is also shifted toward lower wavenumber in the power spectrum obtained from 
the sunspot relative to that from the quiet Sun. For the quiet Sun, convection 
peaks at wavenumber of approximately 2.7 rad Mm"^ that is corresponding to a 
characteristic spatial scale of 2.3 Mm, presumably representing the scale of solar 
granules. For the sunspot, convection peaks at approximately 1.9 rad Mm~^, 
corresponding to a characteristic spatial scale of 3.3 Mm, over 40% larger than 
the typical convection scale estimated for the quiet Sun. This is consistent with 
the result of time - distance analysis of a sunspot observed by Hinode Ca H 
line ( |Zhao, Kosovichev, and Sekii, 2010| ) in that the convection scale is larger 
inside sunspots than in quiet Sun, but different in the convection cell sizes inside 
the sunspot, which was reported as 4 - 5 Mm in that study. The difference in 
sizes obtained from these two studies may be due to the different spectral line 
formation heights. The reason of why the convection cell size appears larger 
inside the sunspot is not clear, but it is possible that magnetic field suppresses 
convection more strongly in smaller scales and leaves larger-scale convection 
more prominent. And also, we cannot completely rule out the possibility of that 
the penumbral filament, which is typically larger than granules, leaks into our 
analysis of convection power inside the sunspot. 
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Figure 6. Comparison of convection power of tiic sunspot and tiic quiet-Sun region. 

4. Discussion 

4.1. Possible Measurement Uncertainties 

Wc have constructed the k~uj power-spectrum diagram using only oscillations in- 
side a sunspot observed by SDO/HMI. Distinct /- and p-mode ridges arc clearly 
seen, and this demonstrates that resonant oscillations exist inside the sunspot, 
primarily inside the sunspot penumbra. It is not clear how much scattered light, 
the light scattered into the dark sunspot umbra from the brighter surrounding 
areas, plays a role in our computed k~uj diagram. It is estimated that in an 
HMI-observed sunspot umbra, the scattered light composes approximately 20% 
of the light intensity (T. Duvall, private communication, 2012). However, the 
power spectrum seen in Figure [2^ is mainly from the sunspot penumbra, and 
this is clear from two experiments that we have performed. One is to mask 
out all oscillations in the umbra, and the other is to enhance the oscillations 
in the umbra to the same amplitude as the oscillations in the penumbra. In 
both experiments, the k-u diagram is not much different from the one shown in 
Figure [5^ except in very low-fc area, where results are already deemed not useful 
due to the selection of a small area and zero-padding. Therefore, the effect of 
scattered light to the computed k-ixj diagram of the sunspot is believed negligible 
in the diagram area that we are interested in. We also attempted to construct 
a k-uj diagram by using only oscillations signals inside the sunspot umbra, but 
were not able to detect distinct modal ridges, probably due to the small size of 
the umbra and its weak oscillation amplitude. 

The uncertainties in determining Doppler velocities inside the sunspot may 
al so have some effects in the k-uj diagram constructed in this studv. Norton et 
al. (|2006p discussed the sensitivitv of the line Fe I 6173A to magnetic field, and 
ICouvidat et a/.l(j2012p examined the uncertainties of deriving Doppler velocities, 
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particularly inside active regions, using different methods. The uncertainties in 
Doppler velocities may cause some uncertainties in power measurements but are 
not expected to play a significant role in shifting the locations of modal ridges 
found in this study. 

4.2. Implications for Sunspot Subsurface Structure 

It is interesting to understand why the p-mode ridges shift toward lower wavcnum-| 
ber (or higher frcaucncv). ISchunker et aZ.I(l2008|) explored the wave motions with 
the presence of inclined magnetic field inside sunspot penumbra, and Felipe, 
Khomenko, and CoUados (|2010p discussed the wave mode conversions in the 
magn etized plasma bv using 3D MHD simulations. More recently. Kitiashvili 
et al. (|201ip showed through MHD simulations that vertical magnetic field 
would introduce some acoustic frequency shifts toward higher frequency, but 
an inclined magnetic field would shift the frequency more substantially. This is 
consistent with our observation of modal-ridgc shifts toward higher frequency 
in the sunspot, where the sunspot penumbra, the region with strong inclined 
magnetic field, contributes most to our computed k-u diagram. It is possible 
that the modal-ridge shift is due to that part of the acoustic waves converts to 
fast magnetoacoustic waves when passing the magnetic field beneath the sunspot, 
and the phase velocity of fast magnetoacoustic wave is determined by: 



where Cs, wa, and are sound speed, Alfven speed, and angle between directions 
of wave-vector and magnetic field, respectively. The observed power spectrum 
is an aggregate from different depths beneath the sunspot, so in principle, it is 
possible to invert the observed power spectrum for the interior magnetic-field 
distribution although this inversion procedure is considerably complicated. 

The relative change in phase velocity (5uph/vph shown in Figure |3J3 may pro- 
vide some information beneath the sunspot. The small values of (Jfph/i'ph, which 
falls in a range of 0.02 — 0.06, together with Equation ([IJ suggests that the 
averaged Scs/cg and v\/c^ inside the sunspot is small, where Scs is the sound- 
speed change inside the sunspot relative to the quiet Sun. If 6cs/cs -C 1 and 
v\/cl <^ 1, from Equation ([!]), we have 



The observed value of Svp^/vph gives a constraint on the value of averaged 
Scs/cs + sin^ ijiv^/c^,), and that is 0.02 < Scs/cs + sin^ ijj{v\/cl) < 0.06. It 
is noted that Scg/cs + sin^ ip{v\/ c^) is averaged over the three-dimensional 
sunspot region weighted with a kernel, which is proportional to the mode energy 
density though the mode energy density inside the sunspot is unknown. Unlike 
siu^ ip{v\/c1) which is always positive, dcg/cs could be positive or negative be- 
neath the sunspot. Since Scs/cg ^ 6j/2j + dT/2T and v\/2c'^ = 1/13, where 7 is 
the adiabatic index and /? the ratio of gas pressure to magnetic pressure. |Jcs/cs| 




(1) 




SOLA: ms.tex; 31 October 2012; 0:43; p. 9 



J. Zhao, D.-Y. Chou 



is smaller than sit? ij^lv^/ c^) near the surface if sin^ -0 is not too small. Thus, 
most of contribution to (5wph/wph could be from siu^ ^{v\/cl). Another possible 
contribution to (5wph/^'phi which is not shown in Equation ([1}, is the change in 
the mode cavity, although it could be smaller than that in Equation ([2|). The 
upper turning point inside the sunspot is deeper than in the quiet Sun, thus its 
contribution to (5vph/vph is positive. 

An interesting phenomenon shown in Figure [3}) is that (Juph/uph has a general 
trend for each n: it first increases and then decreases with k. This could be 
explained by Equation It is generally believed that v\/c?^ is significant only 
near the surface. Since the mode of higher k penetrates less deeply into the solar 
interior, averaged v^/c^ increases with k. On the other hand, the mode of higher 
k propagates less vertically near the surface, and the magnetic field inside the 
sunspot penumbra has a large horizontal component near the surface. This leads 
to that a decrease of "0 with k, or averaged sin^ increases with k. Therefore, 
averaged v\/c1 increases with k while averaged sin^ ijj decreases with k. This 
yields a maximum of sin^ 'ip{v\/c^g) or 5wph/wph at some k. 



4.3. Relevance to Earlier Studies 



Recent power-spectral analysis by ISchunker and BraunI (|201ip showed that the 
modal ridges obtained from the sunspot vicinity, where acoustic halos are de- 
tectable, shift to higher-wavenumber areas relative to the quiet Sun. That shift 
is in a direction opposite to what we find inside the sunspot area in this study. 
It seems that the modal ridges shift toward lower wavenumber in areas where 
acoustic power gets reduced, and shift toward higher wavenumber in areas where 
acoustic power gets enhanced. Moreover, the phase velocity becomes faster in 
the sunspot penumbra possibly implying the presence of faster magnetoacoustic 
waves, and phase velocity becomes slower in the sunspot vicinity possibly im- 
plying a presence of slower magnetoacoustic waves in the area. It is interesting 
to note that in that same study by ISchunker and BraunI (j201ip , the modal- 
ridge lines showed less line-asymmetry than in the quiet-Sun region, again in 
an opposite sense to what is found inside a sunspot in this study. This is also a 
curious phenomenon, but it is not clear whether these opposite line-asymmetry 
trends are related to the different atmospheric heights of sunspot penumbra 
and sunspot vicinity, or related to the faster and slower magnetoacoustic waves 
discussed above. 

The modal-ridge shift observed inside the sunspot may complicate the inter- 
pretation of the measured travel-time shifts or phase shifts in sunspot areas with 
an application of phase- velocity filter or a modal-ridge filter. Because these filters 
choose parameters based on quiet-Sun areas, the power ridges of sunspots may 
not be filtered in the same way as the power ridges of quiet- Sun areas due to that 
the modal ridges of sunspots are systematically located on lower-wavenumber 
(or higher-frequency) side. However, we believe the effect may not be significant, 
as Zhao, Kosovichev, and Sekii (|2010p showed that the measurements with and 



without using phase-speed filters gave results in a good qualitative agreement. 
This observed modal-ridge shift in sunspot regions may not have an effect in 
the ring-diagram derived subsurface properties in active regions, such as given 
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by [Baldner, Bogart, and Basu| (|201ip . because acoustic power of the sunspot 
contributed little to the power spectrum used to do ring-diagram analysis that 
usually analyzes a much larger area than a sunspot. 

5. Conclusion 

By use of 40 hours of continuous SDO/HMI Doppler observations of a stable 
sunspot, we construct a k-uj power-spectrum diagram that shows distinct modal 
ridges with suppressed helioseismic power. More notably, comparing with the k~ 
Lu diagram obtained from a quiet-Sun region, the /-mode ridge in the diagram 
of the sunspot gets more suppressed in power than p-mode ridges, especially at 
high wavenumber. All p-mode ridges shift toward lower-wavenumber (or higher- 
frequency) area for a given frequency (or wavenumber). The phase velocity 
computed from our fitted p-mode ridges shows an increase of 2 — 6% inside 
the sunspot region. Inclined magnetic field in the sunspot penumbra may be 
responsible for the p-modc ridge shifts and the increase of the phase velocity. 
The line-profile asymmetries in the p-mode ridges for the sunspot region are 
more prominent than the asymmetries observed in quiet Sun, and this may 
be relevant to the higher power absorption and suppression ratio in higher 
wavenumbers inside the sunspot. Convection inside the sunspot is also highly 
suppressed, but exhibits a characteristic spatial scale approximately 40% larger 
than the typical granulation scale in the quiet Sun. Realistic MHD simulations of 
magnetoconvection with helioseismic waves included {e.g. [Rempel et al.\ I2009P 
are very useful to understand all of these observed facts inside the sunspot region, 
and we expect that an analysis similar to what is presented in this study but 
using numerical simulation data will shed new light in a better understanding 
of the sunspot subsurface structure. 

Acknowledgements SDO is a NASA mission, and HMI project is supported by NASA 
contract NAS5-02139. We thank Mark Cheung for useful discussions on convection scales 
inside sunspots. We also thank the anonymous referee for useful comments to improve the 
quality of this paper. 

References 

Baldner, C.S., Bogart, R.S., Basu, S.: 2011, J. Phys. Con}. Ser. 271, 012006. 

Basu, S., Antia, H.M., Bogart, R.S.: 2004, Astrophys. J. 610, 1157. 

Braun, D.C., Duvall, T.L. Jr., LaBonte, B.J.: 1987, Astrophys. J. 319, L27. 

Chou, D.-Y., Chen, C.-S., Ou, K.-T., Wang, C.-C: 1992, Astrophys. J. 396, 333. 

Chou, D.-Y., Yang, M.-H., Zhao, H., Liang, Z.-C, Sun, M.-T.: 2009, Astrophys. J. 706, 909. 

Couvidat, S., Rajaguru, S.P., Wachter, R., Sankarasubramanian, K., Schou, J., Scherrer, P.H.: 

2012, Solar Phys. 278, 217. ADS 2 012SoPh..278..217C doi ^l0.1007/sll207-011-9927-y 
Duvall, T.L. Jr., Jefferies, S.M., Harvey, J.W., Osaki, Y., Pomerantz, M. A.: 1993a, Astrophys. 

J. 410, 829. 

Duvall, T.L. Jr., Jefferies, S.M., Harvey, J.W., Pomerantz, M. A.: 1993b, Nature 362, 430. 
FeUpe, T., Khomenko, E., CoUados, M.: 2010, Astrophys. J. 719, 357. 
Georgobiani, D., Stein, R.F., Nordlund, A.: 2003, Astrophys. J. 596, 698. 
Georgobiani, D., Zhao, J., Kosovichcv, A.G., Benson, D., Stein, R.F., Nordlund, A.: 2007, 
Astrophys. J. 657, 1157. 



SOLA: ms.tex; 31 October 2012; 0:43; p. 11 



J. Zhao, D.-Y. Chou 



Gizon, L., Schunker, H., Baldner, C.S., Basu, S., Birch, A.C., Bogart, R.S., Braun, D.C., 

Cameron, R., Duvall, T.L., Hanasoge, S.M., et al: 2009, Space Sci. Rev. 144, 249. 
Haber, D.A., Hindman, B.W. , T oomre, J., Thompson, M.J.: 2004, Sola r Phys., 220, 371. 

ADS I20045oPh..220..37lTn doi |10.1023/B:SOLA-0000031405-52911-08| 

Hill, F.: 1988, Astrophys. J. 333, 996. 

Howard, R.F.: 1990, Solar Phys., 126, 299. ADS I1990SoPh..l26..29"9Fn 

doi:10.1007/BF00153052 
Kitiashvili, I.N., Kosovichcv, A.G., Maiisour, N.N., Wray, A. A.: 2011, Solar Phys., 268, 283. 

ADS:2011SoPh. .268..283K. doi 10.1007/5 11207-010-9679-0^ 
Komm, R., Howe, R., Hill, F., Gonzalez Hernandez, 1., 'Ibner, Tj., Corbard, T.: 2005, Astrophys. 

J. 631, 636. 

Kosovichev, A.G ., Duvall, T.L Jr., Scherrer, RH.: 2 000, Solar Phys. 192, 159. 
ADS 2000SoP h..l92 .159KI doi |10.10 23/A:1005251208431' 

Nagashima, K., Sekii, T., Kosovichev, A.G., Shibahashi, H., Tsuncta, S., Ichimoto, K., Kat- 
sukawa, Y., Lites, B., Nagata, S., Shimizu, T., et al.: 2007, Pub. Astron. Soc. Japan 59, 
S631. 

Nigam, R., Kosovichev, A.G., Scherrer, P.H., Schou, J.: 1998, Astrophys. J. 495, L115. 
Norton, A. A., Graham, J. P., Ulrich, R.K., Schou, J., Tomczyk, S., Liu, Y., Lites, B.W., 

Lopez Aristc, A., Bush, R.I., Socas-Navarro, H., Scher rer, P .H.: 2006, Solar Phys. 239, 

69. ADS:2006SoPh..239...69N, doi: 10.1007/sll207-006-0 279-y[ 
Penn, M.J., LaBonte, B.J.: 1993, Astrophys. J. 415, 383. 

Rempel, M., Schiissler, M., Cameron, R.H., Knolker, M.: 2009, Science 325, 171. 

Rhodes, E.J., Jr., Kosovichev, A.G., Schou, J., Scherrer, P.H., Reiter, J.: 1997, Solar Phys. 

175, 287. ADS:1997SoPh ..l75..287Rl doi|Tai023/A:1004953425123 
Scherrer, P.H., Bogart, R.S., Bush, R.L, Hoeksema, J.T., Kosovichev, A.G., Schou, J., Rosen- 
berg, W., Springer, L.; Tarbell, T.D., Title, A., et al: 1995, Solar Phys. 162, 129. 

ADS:1995SoPh..l62..129S, doi 10.1007/BF00733429 
Scherrer, P.H., Schou, J., Bush, R.L, Kosovichev, A.G., Bogart, R.S., Hoeksema, J.T., 

Liu, Y., Duvall, T.L. Jr ., Zhao, J., Title, A.M., et al: 2012, Solar Phys. 275, 207. 

ADS I2012SoPh..275..207 S doi 10. 1007/sll207-011- 9834-2 
Schou, J., Antia, H.M., Basu, S., Bogart, R.S., Bush, R.L, Chitre, S.M., Christensen-Dalsgaard, 

J., di Mauro, M.P., Dziembowski, W.A., Eff-Darwich, A., et al: 1998, Astrophys. J. 505, 

390. 

Schou, J., Scherrer, P.H., Bush, R.L, Wachter, R., Couvidat, S., Rabello-Soarcs, M.C., Bog- 
art, R.S., Hoeksem a, J.T. , Li u, Y. , Duvall, T.L. Jr., et al: 2012, Solar Phys. 275, 229. 
ADS: 2012SoPh.. 275. .22951 doi [iaT007/sll207-011 - 9842-2 

Schunker, H., Braun, D.C.; 2011, Solar Phys. 2687 349. ADS I2011SoPh.. 268. .34951 
doi:10.1007/sll207-010-9550-3] 

Schunker, H., Braun, TrCT T^E mdse y, C, Gaily, P.S.; 2008, Solar Phys. 251, 341. 
ADS 20085oPh.. 2 51. .34151 doi |10.1007/ sll207-008-9142-7 

Zhao, J., Couvidat, S., Bogart, H.S.," Parchevsky, K.V., Birch, A.C., Duvall, T.L. 
Jr., Beck, J.G., Kosovichev, A.G., Scherrer, P .H.: 2012, Solar Phys. 275 , 375. 
ADS:2012SoPh..275..375Z, doi 10.1007/s ll207-011-9757j^| 

Zhao, J., Kosovichev, A.G., Sekii, T.: 2010, Astrophys. J. 708, 304. 



SOLA: ms.tex; 31 October 2012; 0:43; p. 12 



